Relationship between the Kramers-Kronig relations and negative index of refraction 
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The condition for a negative index of refraction with respect to the vacuum index is established in 
terms of permittivity and permeability susceptibilities. It is found that the imposition of analyticity 
to satisfy the Kramers-Kronig relations is a sufficiently general criterion for a physical negative index. 
The satisfaction of the Kramers-Kronig relations is a manifestation of the principle of causality and 
the predicted frequency region of negative index agrees with the Depine-Lakhtakia condition for the 
phase velocity being anti-directed to the Poynting vector, although the conditions presented here 
do not assume a priori a negative solution branch for n. 



Negative index optics is an attractive field because it 
is driven by a wealth of exotic applications such as the 
creation of a perfect perfect lens \V\, cloaking i2] and 
through- wall vision |3|]. While much success has been 
made of actualizing a negative refractive index in 'meta- 
materials', there are still some theoretical loose ends in 
the understanding of the physical conditions in which 
negative refraction can occur, quite apart from consid- 
erations of efficiency losses, material choice and exper- 
imental proof of principle. Experiments by Shelby Q 
and Smith [&\ have shown negative refraction in the mi- 
crowave region of the electromagnetic spectrum using 
meta-materials, while metal/dielectric hybrid rod struc- 
tures allowed Yao [6|| et al. to demonstrate negative re- 
fraction in the visible spectral region. 
The experimental proof of negative refraction lies in the 
arena of geometric ray bending in slabs and prisms as 
well as being derivable from the transmittance and re- 
flectance data. The original criteria which were postu- 
lated for a material to possess a negative index of re- 
fraction (n=Re[N]) were given by a paper by Veselago 
Q in 1968, wherein he established that, in order to have 
a negative index, the real parts of e and /^ need to be 
simultaneously negative. In this scenario, the magnetic 
induction, electrostatic displacement and the wave- vector 
form a left-handed triad and the material possesses a neg- 
ative index in the frequency region wherein the Veselago 
condition is satisfied. However, it has emerged that the 
Veselago conditions are not sufficiently general to cap- 
ture all eventualities of negative n, because they are only 
valid for a purely real permittivity and permeability. An 
important paper by Depine and Lakhtakia [8j showed 
that, having assumed the negative branch of the solu- 
tion for n, the condition of having the phase velocity 
anti-directed to the group velocity (and Poynting vec- 
tor) gives a more general criterion for negative refraction. 
This was a generalization of the concept of negative re- 
fraction to negative phase velocity. Nevertheless, there 
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exists disagreement in the literature as exemplified by the 
theoretical presentation of a case by Valanju et al. [9] of 
positively refracting waves whose group and phase veloc- 
ities were non-collinear in an inhomogeneously dispersed 
wave packet in a negative index material. 
McCall [10| considered a space-time covariant formalism 
of the calculation of phase velocityin negative index ma- 
terials. Mackay and Lakhtakia [ll| recently considered 
the effect of bianisotropy in determining the directional- 
ity dependence of negative index, phase and group ve- 
locity conditions in terms of the complex wave- vector for 
pseudo-chiral omega materials. In this work, we consider 
scalar responses, but we examine the model independent 
criteria for negative index in terms of the response func- 
tions. Stockman [12| considered the square of the com- 
plex refractive index N^ and applied the Kramers-Kronig 
(KK) relations (if e and ^ are causal response functions, 
N^ is analytic and also obeys the KK relations) to relate 
the phase and group velocities to the dissipation in the 
system. This approach worked with complex N^ instead 
of the real n and k. Since the latter two are the optical 
constants which determine the phase velocity and group 
velocity, and plane wave attenuation, it is the purpose 
of this letter to work with these and impose the KK re- 
lations at this level. Kinsler ^3|] finds dealing with the 
refractive index under the definition n—y''Re{N^) more 
convenient because one does not need to make correc- 
tions for spatial oscillations of the counter propagating 
electric field at interfaces. Nevertheless, there is still a 
sign choice implicit in this definition and he appeals to 
Dephine and Lakhtakia's work J8|] for the criterion for 
negative n which e and fi must satisfy. Peiponen [14l | 
showed that, having solved for the refractive index us- 
ing the Veselago conditions, the real and complex parts 
of the refractive index satisfied the Kramers-Kronig re- 
lations. The purpose of this letter is to show that the 
satisfaction of the Kramers-Kronig relations should be 
the primary criterion for determining the correct sign of 
n. 

The above considerations suggest a need of a more ro- 
bustly general way of looking at the conditions for nega- 



tive n with respect to the vacuum refractive index. In this 
paper, we wish to estabhsh the most general criteria for 
negative n (Re[N]), by examining the solution branches 
for N^=e/i, and choosing the solution in each frequency 
region which maintains analyticity (differentiability along 
the real axis of n). The choice of the solution such that 
the functions ti{uj) and k{uj) (the extinction) are analytic 
implies that these functions, being the real and complex 
parts of a function N=n+iK, must satisfy the Kramers- 
Kronig relations [la]. This is an indirect manifestation 
of causality at the level of electromagnetic plane waves - 
there can be no dispersion without absorption. 
Let us examine all possible solutions to N^ = (e/i), where 
the permittivity and permeability are given by e and /i. 
Writing e=er+i Ci and using a similar expression for /i, 
we arrive at the following quartic : 



CrfJ'i ~r CifJjr 



{erfJ-r - eitJ-i)n'^ = 0. (1) 



A similar equation exists for k. There are four solutions 
for n, two of which are purely imaginary and two of which 
are real and these two are labeled n_|_ and n_. The two 
real solutions are written as follows : 



n± = ± 



V2 



(2) 



These solution branches are plotted in Figure [U and it 
can be seen that there are two separatrices. For this 
plot, we have used one model whereby a Lorentzian res- 
onance is used for both the permittivity and permeabil- 
ity (Lorentz-Lorentz), while the second model assumes a 
Drude plasmon model for the permittivity response and a 
Lorentzian for the permeability (Drude-Lorentz model). 
The models are described in the following way : 
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where the parameters used for the model calculation are 
chosen to be aJom=0.5, woe^l-lS, Wpe^l.O, Wpm^O.S, 
7e=0.3, 7„i=0.11, 600=0.4 and a frequency scale of a;o=l- 
The locations of the separatrix in Fig [T] are two points 
in frequency space at which the curves for nj- are no 
longer differentiable (analytic). We can see this clearly 
by calculating dn-\.ldio and this is written as follows : 
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FIG. 1: (Color Online) Plot of the family of solutions for n 
for which a Drude model is used for e and a Lorentz model is 
used for /i. 



where /i indicates, for example, differentiation with re- 
spect to w. As can be seen from Figure [21 the pairs of so- 
lutions (n_,K_) and (n+,K+) do not satisfy the Kramers- 
Kronig relations. This is because, as can be seen from 
the equation for dn^/duj above , the solution branch 
n+ departs from analyticity (by possessing a divergent 
derivative) at the frequency which satisfies the following 
equation : 



e^Mr - ^^^J■l + |e||Ai 



= 0. 



(6) 



Mathematically, n_|_(a;) and n_(a;) are not analytic be- 
cause of presence of the absolute value function. 

The solutions to Equation |6] define the locations of the 
separatrix in Fig [1] The departure of analyticity here 
means that, without the correct sign choice, the group re- 
fractive index would be divergent and un-physical. Not- 



ing that, |e| 
read : 



^ 



we can simplify Equation [6] to 



Mi 



(7) 



This equation admits at least two solutions ui± for the 
Lorentz-Lorentz and Drudc-Lorentz models used here, 
and the number of solutions in general depends on the 
model for e and /i used. The condition above in Equa- 
tion [7] is actually identical to that which was suggested 
by Ruppin [16] in a book chapter jl.7,j and previous article 
[18| . The condition is equivalent to the Depine-Lakhtakia 
condition [8[ except for a sign ambiguity. If we choose the 
solution for n such that 

n+(a;), < a; < w_ 
n{uj) = |n_(w), a;_ < oj < a;+ 
n+{uj), w+ < w < oo 
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FIG. 2: (Color Online) Plot of the n+ solution branch, 
together with its Kramers-Kronig transformation for the 
Lorentz-Lorentz model (a) and the Drude-Lorentz model (b). 
Note that the Kramers-Kronig transformed n does not match 



which now satisfies the Kramers-Kronig relation [19| : 



k{uj) 



{n{uj ) — l)w 



du 



which can be seen from Figure |4l This scheme for choos- 
ing the sign of the refractive index requires the frequency 
range to be calculated by Equation [71 and the phase an- 
gle of N^ (=(n-|-iK)^) lies in the within the interval [-7r,0] 
in the frequency range w_ < w < a;+, while it lies in [0,7r] 
everywhere else. For the Lorentz-Lorentz case computed 
here give w_=1.21 and a;+=1.43 in units of wq. 

These are the frequency regions for which we must 
switch to the negative branch n_, in order to main- 
tain differentiability in n(a;) . This now defines a set 
of (n(u;),K(w)) which satisfies the Kramers-Kronig (KK) 
relations, as plotted in Fig 2] (a). This requirement of an- 
alyticity to satisfy the KK relations is a manifestation of 
the relationship between the real and complex parts of an 
analytic function and is a manifestation of the principle of 
causality. It is shown in Figure [3] that the Veselago condi- 



FIG. 3: (Color Online) (a) Plot of the refractive index un- 
der the Veselago condition (choosing the n<0 solution branch 
(n_) when e and ^ are simultaneously negative) for the 
Lorentz-Lorentz model, (b) a comparison between that Vese- 
lago refractive index solution and the solution determined by 
the Depine-Lakhtakia condition. We see that that Veselago 
criterion is not sufficiently general, because it does not sustain 
the validity of the Kramers-Kronig relations. 



tions are not sufhciently general to capture the frequency 
interval over which the refractive index is negative. If we 
were to choose the n_ only in the region where Re(e) and 
Re(/i) were both less than zero, we would violate analyt- 
icity. 

It can be seen from Figure H] that the predicted fre- 
quency region of negative n agrees with those given by the 
Depine-Lakhtakia condition 0], which was established 
from considering that the phase velocity vector should 
be anti-directed to the Foynting vector. This comes as a 
corollary to our Kramers-Kronig criteria here, most sim- 
ply encapsulated in Equation [T] 

We have found that the principle of causality alone, as 
expressed in the Kramers-Kronig relations, is sufficient 
to establish whether a material possesses a negative in- 
dex of refraction with respect to vacuum, and if there is 
a negative index, the frequency range for which the re- 
fractive index is negative are given. In frequency space, 
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causality implies that the refractive index is analytic in 
the upper-half complex plane, unless the responses of the 
system are active, wherein the condition of analyticity 
is sustained by choosing the integration contour for the 
Kramers-Kronig relations such that poles lie on a line be- 
low this contour, as described by Skaar [20] ■ The analyt- 
icity condition is consistent with those of the Dephine- 
Laktakia analysis, but we don't not need to assume a 
■priori that n <0. In this sense, the choice of the negative 
branch based on sustaining analyticity of n(w) and k[lS) 
(and hence to satisfy of the Kramers-Kronig relations) is 
the most general criterion. 
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FIG. 4: (Color Online) (a) a comparison between the 
Depine-Lakhtakia condition and the analytic condition for the 
Lorentz-Lorentz model, (b) The analytic condition imposed 
on the solution branch restores the Kramers-Kronig relations 
, and hence causality. 
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